Flux transitions in a superconducting ring 
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We consider a superconducting ring embedded in a magnetic field. The magnetic field is initially 
too high to support a superconducting state, but is ideally slowly decreased and then reversed. 
During this process, the number of fluxoids trapped by the ring exhibits a complex pattern of 
transitions, which depends qualitatively on the ratios of the geometric dimensions of the ring to 
the coherence length. We evaluate the current around the ring and find rough agreement with 
experiments. However, we also predict a new feature that has not shown up in experiments: in the 
vicinity of a critical field, giant jumps are expected. 
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I. INTRODUCTION 

It has been long known that a superconducting loop 
embedded in a magnetic field exhibits phenomena of fun- 
damental nature, such as persistent currents and flux 
quantization.EH3 The concept of flux quantization has 
been crucial for the analysis of a wide range of topics 
of basictH and practicaltl importance. The ability to 
manufacture loops with widths of the order of the co- 
herence length- brought renewed interest in this system; 
experimental]!] and theoreticalEM studies are being per- 
formed until present. 

The purpose of the present article is the theoretical 
study of the following situation. A ring of supercon- 
ducting material is embedded in a homogeneous magnetic 
field. Initially the field is too high to support a supercon- 
ducting state, but it is slowly lowered to zero while the 
temperature is kept constant; the field is then reversed 
and increased until superconductivity is destroyed. The 
most prominent feature of these experiments comes from 
the thermodynamic need for a winding number (num- 
ber p£ fluxoids) which is not too far from the applied 
flux,E9 together with the constraint that this number has 
to be integer. These coupled requirements usually lead to 
jumps in the winding number and, accordingly, in quan- 
tities such as the current around the ring. 

Experiments like this were performed long .agn.. but 
modern versions of it are still giving new results. EiE3ll3'0 
Situations similar to that considered here have been re- 
ported in Refs. [jj^,[l6|. In the limit of a thin ring, this 
situation provides a neat exanialp-|for the study of state 
selection in phase transitiopsJl3'L3 Similar .experiments 
for weakly connected rings,E2l and for disksE3'ErE3 have 
also been widely studied. 

In the limit of infinite radius, a ring becomes a slab 
with parallel plane boundaries. We shall suggest in 
Sec. [n] that the most outstanding feature obtained in 
our results can be related to the condition for the for- 
mati on , of .vortices in a slab parallel to the magnetic 
field £aEJE§E£l If the coherence length is larger than 0.55 
times the width of the slab, then vortices cannot be 
present in the slab, whereas for smaller coherence lengths 



they can. This change of regime is accompanied by a 
change of behavior of the temperature-field dependence 
along the normal-superconducting (N-S) phase bound- 
ary. The second derivative of the transition temperature 
with respect to the field is discontinuous at the criti- 
cal value of the coherence-, length. These results were 
confirmed exp er iment ally.EH For more recent studies on 
slabs, see e.g. Refs. p8|]29| , |30| . 

We consider a ring with perfect axial symmetry and an- 
alyze it by means of the Ginzburg-Landau model. Since 
the magnetic field varies with time, it is natural to use 
the time-dependent Ginzburg-Landau model (TDGLX. 
This approach has been followed in recent calculations !El 
While TDGL is appropriate for the study of a given 
experiment, the number of parameters it introduces is 
rather large, depending on the material, the tempera- 
ture, and the rate at which the magnetic field is swept. 
In this article we have chosen to consider the theoretical 
limit in which the field varies ideally slowly, so that the 
time independent formalism can be used. 

We assume that the ring is very thin compared with the 
magnetic penetration length. Under these assumptions 
the free energy is (up to unimportant constants) 



G= I <M I ,-,/,- M (-|^| 2 + i|^|4) + |(i#V-&r/#0)Vf 

where r and 9 are cylindrical coordinates, ip is a normal- 
ized order parameter, i?i and R are the inner and outer 
radii of the ring, /i = (i?/£) 2 is a function of the tem- 
perature (£ is the coherence length) and b — itR 2 H c /& , 
where H c is the magnetic field and $q the quantum of 
flux. 



II. THE ONSET OF SUPERCONDUCTIVITY 

In order to gain some intuition about this problem, 
we first consider the Little-Parks regime, at the onset of 
superconductivity. In this regime \ip\ <C 1, so that the 
Ginzburg-Landau equations are liueapzed and the case 
of the ring can be solved exactly.OOo The purpose of 
the present section is to point out that there exists a 
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critical point in the N-S boundary. As will be seen in 
the following sections, it appears that this critical point 
plays a major role in the "superheating" stability of the 
superconducting states. 

Figure [l] shows the normalized field b at which a ring 
with i?i = 0.8R becomes superconducting, as a function 
of the winding number to, for fi = 92 and for fi = 94. 
Only integer values of to are meaningful, and only the 
value of to with the highest b is physically realized, but 
these curves are useful for explanatory purposes. We 
found that for /i < 92 these onset curves have only 
one maximum, whereas for /.i > 94 they have two. At 
/i ~ 93, to ~ 68, b ~ 87 there is a critical point, where 
the first and the second derivative of the onset field (or 
of the free energy) with respect to to, both vanish. As 
a consequence, near this critical point, the free energy 
depends very weakly on the winding number to. 

The influence of this critical point is felt quite far from 
its immediate vicinity. The dashed line in Fig. |l| is a 
line of constant (negative) free energy for /j, = 110. We 
see in its upper part that for a range of about 8 winding 
numbers the slope of b with respect to to, expressed in 
"natural" units, is less than a tenth. By the relationship 
(db/dm)a = — (dG/dm)i,/ (dG/db) m , the smallness of 
(db/dm)c suggests that the free energy also has a very 
weak dependence on the winding number. 

In reality, only integer values of to and their linear 
combinations exist. The typical scenario for flux tran- 
sitions is as follows. A given field b is present and the 
order parameter has a certain winding number to. Lin- 
ear combinations with other winding numbers act as en- 
ergy barriers and the sample is "trapped" at m. When 
b becames sufficiently small, the energy barrier towards 
smaller winding numbers is levelled off, and the state de- 
cays. If (dG/dm)b is large, the order parameter will be 
"stopped" by the following barrier and become trapped 
at to — 1. However, a small slope means that to has to 
change by a big amount for a given change in G, leading 
to the possibility of transitions that change to by several 
units. 

The critical point found here can be associated to two 
phenomena which have been discussed in the literature. 
The first is the condition for the formation of vortices in a 
slab, quoted in Sec. |. Using the value 0.2i? as the thick- 
ness of the slab, the onset for the appearance of vortices 
is located at (i ~ b ~ 80, not far from the critical point in 
our case. We can therefore expect that for n above the 
critical point (low temperatures) there will be vortices in 
the ring, whereas for small fi the winding number will be 
the same for the inner or for the outer boundary of the 
ring. 

An additional feature noted in Ref. ^4] is the crossover 
from a regime in which at the N-S boundary b is approx- 
imately proportional to fi 1 / 2 to a regime where b is ap- 
proximately proportional to /i. In the case of rings, this 
crossover applies to the background line of this boundary, 
which-has the Little-Parks oscillations superimposed on 
itBJ'L3 This is actually a dimensional crossover: at low 
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FIG. 1: The continuous lines mark the onset of superconduc- 
tivity in the flux versus number of fluxoids plane. The ring 
is superconducting in the area "inside" the appropriate line. 
Only integer values of m are meaningful. The ratio between 
the inner and the outer radius is 0.8. The lower line is for 
(j, = 92 and the upper line, for \i = 94. The dashed line is 
a line of constant energy for [i — 110. The continuous lines 
were obtained by an exact calculation; the evaluation of the 
dashed line approximates the radial dependence of the order 
parameter by a sum of three orthogonal polynomials. 

fields the sample behaves as a thin slab and b = 5(3/z) 1 / 2 ; 
at high fields it behaves as half a plane and b = 0.847/1. 
Again, the crossover occurs at /x ~ 100. For the case of 
slabs the passage from one regime to the other involves 
discontinuities of some quantities, but for the case of rings 
these discontinuities are not relevant, since except for ac- 
cidental situations the critical point does not occur at a 
physical (i.e. integer) value of to. 

Since we believe that the main features of the situation 
we study are governed by the critical point described in 
this section, we have evaluated numerically the position 
of this point as a function of the ratio R-JR. Figure || 
shows the values of fi and b at which the critical point 
occurs. Since to a rough approximation these values are 
inversely proportional to the square of the width of the 
ring, the normalized values n(l — R\/R) 2 and b(l — RjR) 2 
appear in the figure. In the following sections we will 
always take R\ = 0.8R, which is a representative experi- 
mental value. 



III. MATHEMATICAL METHOD 

We return now to the general situation in which the 
order parameter -0 is not necessarily small. We write tp as 
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FIG. 2: Normalized values of fi and b at which the critical 
point occurs, as functions of the ratio between the inner and 
the outer radius of the ring. The upper curve is — RJR) 2 
and the lower curve is 6(1 — Ri/R) 2 . In the calculation of these 
curves we have approximated the radial dependence of the 
order parameter by a sum of three orthogonal polynomials. 
For Ri — * R, the critical values of /i and 6 diverge and our 
method cannot be used; in this limit, r+he ring becomes a slab 
and m expect fi(l - Ri/R) 2 -> 3.393 and 6(1 - Ri/R) 2 -» 
3.24.E3 



a function of the cylindric coordinates r and 6 and assume 
that the ring is so thin that there is no z-dependence. We 
expand 



(r)e~ 



l P n (r) , 



(2a) 



(2b) 



where Eq. (ga|) is a Fourier expansion and Eq. (Eq) is an 
expansion into orthogonal polynomials, with P^\R{) = 
Pn(R) = and such that P n (r) vanishes n times in the 
range Ri < r < R. We then keep up to four terms in 
Eq. (^aj) : two "leading" terms with winding numbers too 
and toi , and two "satellite" terms, with winding numbers 
2too — toi and 2mi — too. In Eq. (2b) we keep the terms 
n = 0, 1, 2 for the leading terms, and only the term n = 
for the satellite terms. The coefficients c mn which are 
not declared to vanish, are found by minimizing the free 
energy. 

The physical meaning of a single winding number in 
the expansion (^a|) is that has axial symmetry. If 
there are two leading winding numbers Too and Toi and 
two satellites, as we shall consider as our most general 
situation, and if the weights of the two leading winding 
numbers are similar, then the equilibrium order param- 
eter exhibits a ring of |too — toi equidistant vortices in 
the superconducting sample. If there arc several winding 
numbers, but the weight of one of them is much larger 
than those of the others, then \ip\ does not have axial 
symmetry, but there are no vortices in the sample. For 
situations in which there are vortices in the sample that 



FIG. 3: Comparison of approximated and exact values of the 
order parameter for a typical situation. The difference be- 
tween both curves does not exceed 0.001. 



are not equidistant from the z-axis, our approximation 
breaks down. 

According to our results, the approximations used 
seem to be reasonable as long as /i is not much higher 
than that of the critical point discussed in Sec. ||. The 
usual situation is that the weight of the leading terms 
is much greater than that of the satellites. In Fig. |^ we 
compare the radial dependence of an axially symmetric 
order parameter obtained by the present approximations 
to that from an exact calculation. 

Our procedure is as follows. For fixed temperature 
(fixed /i), we first evaluate the magnetic field for the on- 
set of superconductivity. At the onset, the Ginzburg- 
Landau equations are linear, so that the order parameter 
generically has axial symmetry and the sum (|2a|) con- 
tains only one to, which we denote by too. The field is 
then gradually lowered and for every field we evaluate 
the order parameter lZ mo (r)e~ m ° lS . We then check the 
stability of this order parameter against small variations. 

Variations explore subspaces with larger dimensions 
than that of the equilibrium order parameter. When 
looking for decays from too, variations also contain the 
winding numbers too — I and too 4- £, where £ is swept 
over at least the range 1 < £ < n/5. It is important to 
include both too — £ and too 4- £, since due to nonlinearity 
of the GL model these harmonics are coupled. For each 
variation we evaluate the Hessian of the variation of the 
free energy and obtain its eigenvalues. If all the eigenval- 
ues are positive, it means that lZ mo (r)e~ m ° ie is stable, 
and we proceed to decrease the field by a small step. 

If any of the eigenvalues is negative, it means that 
Kmo( r )e~ m °'" is unstable with respect to the variation 
represented by the corresponding eigenstate. We denote 
by toi = too ± £ the value of the winding number that 
has the highest weight in the variation against which too 
was unstable. 

During the decay of TZ mo (r)e~ m ° ie , the order parame- 
ter changes quickly. Therefore, in order to follow the or- 
der parameter or the currents during the decay, we would 
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have to use the TDGL model. Using a gauge in which 
the electrochemical potential vanishes, the equations of 
this model take the form 

d^/dt = -D^G, (3a) 



dA/dt = -D A S A G , (3b) 

where A is the vector potential, 5^ and i5a denote vari- 
ations over the complex conjugate of tp and over A, and 
and D A are positive constants. From here, it follows 
that G will always decrease during the decay and the or- 
der parameter will reach equilibrium when G reaches a 
minimum. (The applied field is assumed to remain con- 
stant during the decay.) 

Therefore, in order to obtain the order parameter after 
the decay, we perform a minimization of the free energy 
over states of the form (^), where the leading winding 
numbers are Too and m\. The initial state from which the 
minimization flows is a combination of the former state 
with to = mo and of the eigenstate that had the nega- 
tive eigenvalue. The minimization may lead to a state in 
which only one winding number has a significant weight; 
in this case we recover the previous situation and we just 
have to update the value of mo. The other possibility is 
to be left with a combination of winding numbers (which 
means that \ip\ is not axially symmetric). It seems plau- 
sible to assume that in the generic case the minimization 
leads to the same minimum, and with the same leading 
harmonics that would be obtained by TDGL; we intend 
to investigate the range of validity of this assumption in 
a future study. 

In any numerical procedure, arbitrary decisions have 
to be taken concerning parameters such as cutoffs and 
convergence factors. Usually these parameters are fixed 
empirically in order to optimize convergence properties. 
In our case, one has to choose the initial weights given 
to the original state and to the deviation that renders 
it unstable in their initial combination, the minimization 
method, and the minimum weight below which a winding 
number is considered "not significant" . We have used 
several criteria and become convinced that, in most cases, 
our results do not depend on the particular choices. 

If the order parameter is of the form (^a|) with lead- 
ing winding numbers too and mi, we check its stabil- 
ity along similar lines to those described for the axial- 
lly symmetric case. This time the attempted variations 
contain combinations of winding numbers £, 2too — £ 
and 2toi — £, where £ is swept over at least the range 
[max(TOo, toi) — /i/5] < £ < [max(mo, toi) — 1]. 

IV. RESULTS 

For each field, after the order parameter is known, the 
free energy is evaluated using Eq. (Q) , the normalized cur- 
rent density is evaluated as Re[^(ii?V — br/R8)i(j] and 
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FIG. 4: Normalized current as a function of the normalized 
field for fi — 0.475 and for fi = 0.95. Decays are repre- 
sented by the vertical lines at the left side of each undula- 
tion. For visibility, the curve for fi — 0.95 has been raised 
by 0.06. The physical current per unit height, in electrostatic 
units, is obtained by multiplying the normalized current by 
c$oM/2(27TKii) 2 , where c and k are the speed of light and the 
Ginzburg-Landau parameter. 

then the current (per unit height) is obtained by integra- 
tion over r of the tangential component. 



A. Near T c 

Figure || shows the current around the ring, as a func- 
tion of the field, for low values of fi. As expected, 
there are discontinuities in the current, nearly period- 
ically spaced. These discontinuities correspond to the 
passage between consecutive winding numbers. 

For /i = 0.475 and close to b = ±4.35, the ring is 
in the normal state. This reentrant behavior is due 
to the Little-Parks oscillations. For a ring with Ri = 
0.8R, reentrant superconductivity is possible in the range 
0.02 < n < 12 -4. For ^ < 0.02, only the Meissner 
state (to = 0) exists; in the range 0.02 < [i < 12.4 
there are "windows" for which superconductivity is in- 
terrupted when passing from one to to the next. For low 
values of /j, these windows are wide, but, as [L increases, 
the background slope of the N-S boundary in the b — fi 
plane increases and these windows become narrower. For 
/i > 12.4 the N-S boundary is monotonia 



B. Short coherence length 

Figures |S| and |6| show our results in the range 20 < 
fi < 120. For [i > 120, our approximation of Eq. ( |2a| ) by 
four uniformly spaced harmonics breaks down. This was 
expected since, if the distance between the boundaries is 
considerably larger than the coherence length, then there 
is room for more than one ring of vortices. 

Obviously, if the magnetic field is reversed, the current 
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FIG. 5: Normalized current as a function of the normalized 
field for a wide range of values of fi. The value of n is written 
next to each curve. Shifts of 0.1 have been inserted between 
consecutive lines. The right end of this graph is shown en- 
larged in Fig. ^. 
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FIG. 6: Free energy as a function of the normalized field for 
the same values of n as in Fig. || Shifts of —2 have been 
inserted between consecutive lines. The central part of fj, = 
120 has been chopped. The right end of this graph is shown 
enlarged in Fig. 



in the sample for the equilibrium (absolute minimum) 
configuration is reversed too. For small values of /x hys- 
teresis is a minor effect, and indeed we note in Fig. |] 
that the current for /i = 0.475 is nearly an odd function 
of the field. On the other hand, for the values considered 
in Fig. |^, hysteresis is dominant and the equilibrium val- 
ues have practically no importance. At first glance the 
curves in Fig. |^ look like even functions and we observe 
that qualitatively different regions appear in nearly sym- 
metric pairs. In the following, we shall review these re- 
gions, by decreasing values of \b\. Following the line of 
reasoning of Sec. ||, we may argue that the approximate 
evenness is due to the invariance of \(dG/dm)i,\ when 
both b and m change sign. 



FIG. 7: Normalized current near the onset of superconductiv- 
ity. The lowest line is for fi = 110 and sits at its true position. 
The other lines are for /i = 120 and fi = 130. For each in- 
crement of \i by 10, the position of the line is raised by 0.015 
and is also shifted by 7 to the left. 
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FIG. 8: Free energy near the onset of superconductivity. The 
steepest line is for fi = 110 and sits at its true position. The 
other lines are for /i = 120 and [i = 130. For each increment 
of fjj by 10, the position of the line is shifted by 7 to the left. 



1. Near the onset 

Figures [j] and |^ show our results for high values of b, 
with /i = 110, 120 or 130. [Near the onset of supercon- 
ductivity, the approximation of the expansion (2a) by a 
limited number of terms remains valid beyond /i = 120.] 
The region at the right consists of steps where a single 
value of m is left in (Eat) and, as b is lowered, m decreases 
by 1 between consecutive steps. There are five such steps 
for [i = 120 and 130, and four steps for /i = 110. The 
current always drops when the state decays to a lower m. 



2. Breaking of axial symmetry 

This is roughly the region 87 < b < 94 for [i = 110, 
81 + 7 < b < 94 + 7 for /i = 120 and includes most of 
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the line for fx = 130 in Fig. |t]. If we denote by (mo) a 
state with a single harmonic to — m Q and by (mo, mi) 
a state with two leading winding numbers mo and mi, 
this region extends from (84, 70) to (77,65) for ll = 110, 
from (91, 74) to (78, 66) for fx = 120 and from (98, 79) to 
an uncalculated state for fj, = 130. 

We may divide this region into a subregion of "incip- 
ient symmetry breaking" and a subregion of "developed 
symmetry breaking" . The subregion of incipient sym- 
metry breaking consists of the sequence of double steps 



(m ') (m 

(i+i) 



Here to l 

/ 0) 

(m J ,m\') 



(1) m«) 
Ai+l) 



(m ') -> (m 
,W 



(2) (2) 
' m l 



1 > m\ , so that the transitions 
(ttiq^ 1 ') are discontinuous; on the other 



hand, the transitions (m^) — > (nig 1 ' ,m\ l> ) are continu- 
ous bifurcations. In the subregion of developed symmetry 
breaking, all the states are combinations with two lead- 
ing winding numbers. For /i = 130 there is only one step 
with incipient symmetry breaking and for fi = 120 there 
are two steps; for fi = 110, incipient symmetry breaking 
extends over the entire region. Some segments are very 
short, and cannot be seen in the graph. 

The subregion of developed symmetry breaking can be 
subdivided further: in its high-field part the transitions 
(to^to^) — y (mg 1 ,m[' + ) have either m^ 1 ' = ml'' 

or mj = )Tij , for lower fields there may be "cas- 
cades" . By a cascade we mean that, after a state decays, 
the new state is already unstable with respect to a third 
state, and decays at the same field at which it appeared. 
By means of cascades, too and mi can both change si- 
multaneously. 

For axially symmetric states, the slope of the current 
in Fig. [j] is negative. On the other hand, for combina- 
tions with two leading winding numbers, the current is a 
markedly convex function, and can reach large positive 
slopes. For decays that follow after these large positive 
slopes, the current usually rises. 



3. Critical region 

This is the region with large discontinuities near |6| ~ 
70. For (j, = 120 and fi = 110 the region of asymmetry 
ends with a long cascade, the last stage of which is a 
discontinuous decay into a symmetric state. For /1 < 100 
the region of asymmetry is absent. For 11 = 100 the 
critical region begins with a decay in which m overshoots 
and then goes up. For b < 0, the critical region in which 
to jumps by several units appears without previous signs. 

It may be surprising that although the samples are 
sufficiently wide to have vortices on them, the order pa- 
rameter is axially symmetric in this region. An intuitive 
explanation might be that the field is not sufficiently large 
and therefore a ring of vortices is not thermodynamically 
favorable. 

Table Q is a record of the states through which the 
critical region passes. 

One notices that for b > the number of fluxoids that 



TABLE I: Fluxoid numbers occurring in the critical region. 
Negative numbers are for the region in the negative field. 



/' 



states 



120 (61) -> (53) -> (48) 

120 (-39) -> (-41) -> (-45) ->■ 

(-51) -> (-59) -> (-69) -> (-79) 

110 (66) -> (57) -> (51) -> (48) 

100 (74) -> (64) -> (56) -> (51) 

100 (-43) -> (-45) -> (-49) -> (-55) -> (-63) -> (-72) 

95 (75) -> (65) -> (58) -> (53) -> (51) 

90 (72) -> (64) -> (57) -> (53) 

80 (65) -> (62) -> (59) -> (57) 

80 (-51) -> (-53) -v (-56) -> (-60) -> (-64) 



go through the ring during the critical region is largest 
for fi near its critical value. The largest signal in the 
current-field graph, as well as the number of fluxoids 
that go through for b < 0, appear to be at higher values 
of Li. 



4- Subcritical region 

This is the broadest region in Figs. || and ||. As the 
field decreases through the critical region, jumps in the 
winding number decrease from giant to moderate, until 
we finally reach a region in which the winding number 
changes by just one in each decay. For fi < 60 the two 
previous regions are absent and there is no distinction 
between this and the first region. 

For 20 < /i < 60 and b > 0, the average range Ab 
between consecutive jumps in this region is 1.30; for b < 
0, Ab = 1.17. This difference can be understood, since for 
b > a delay in the decays is produced by "supercooling" 
and in the negative region it is necessary to catch up with 
this delay. For 80 < Li < 120, Ab = 1.24 for positive fields 
and Ab = 1.23 for b < 0, indicating that the presence of 
the critical region helps to release most of this delay. Note 
that Ab is not close to unity due to our choice of the outer 
radius as the unit of length. Had we chosen the average 
radius 0.9i? as the unit of length (as in Ref. O), we would 
have had to divide the values of b by 1/0. 9 7 « 1.235. 



5. Low fields 

This is a short region centered at slightly negative 
fields. In this region the "satellite" harmonic 2toi — mo 
is not necessarily small. In this region the results we ob- 
tain are sensitive to the minimization strategy, and the 
results in Figs. || and || are less reliable here than in the 
other regions. 
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FIG. 9: Comparison between our results and those measured 
in Ref. [w[ The experimental data are represented by the 
dense set of dots and the theory by the line above. The ex- 
perimental value of b was obtained by multiplying the applied 
field by irR 2 /<f> . Also, since in our calculations the field was 
swept downwards and in the experiment upwards, the sign of 
the experimental field was switched. The experimental value 
of the normalized current was obtained by multiplying the 
field induced by the supercurrent at the center of the ring by 
27TK 2 (7?i + R)R 2 /$ofJ-t. We took k ~ 0.5 and estimated the 
field induced at the center of the ring as 2.5 times its average 
value on the active square of the Hall cross. The calculated 
curve is for n = 30. 



V. COMPARISON WITH EXPERIMENTS 

To our knowledge, the experiments closest to the situ- 
ation discussed in this article are Ref. [ll] and Refs. |l3|Jl4|. 
Reference [ll] deals with only three fiuxoids, and from the 
present point of view may be considered featureless. (See 
however Ref. p6[ ) The results of Ref. [l4| are similar to 
those of Ref. |13|, but in this case the susceptibility goes to 
zero at the onset of superconductivity, as in our results. 
Ref. |lj also confirms our result that, if the width of the 
ring is appreciably smaller than the coherence length, all 
the transitions involve the passage of just one fluxoid. 

Reference [51] reports on two rings. In both cases the es- 
timated mean radius is 2.16/xm and the coherence length, 
180nm. For their wide ring, R — iij = 630nm. Using 
graph ||, these values result in /i = (i?/£) 2 ~ 3jU cr i t , 
far beyond the region we have studied. One observation 
can nevertheless be adventured: the critical value of the 
field is 6 cr it ~ 56, which for this geometry corresponds to 
about GOG. The "catastrophic behavior" in their Fig. 5 
begins not far from there. 

Their narrow ring has a width R — R\ = 316nm, which 
now implies [i — (i?/£) 2 ~ 0.9/i cr it. Translated to our 
case, this corresponds to fi ~ 80. Inspection of their 
Fig. 2 and our Fig. || suggests that fj, = 60 would give 
a better comparison. The difference between /i = 80 
and \x = 60 could be due to a 15% error in any of the 
reported lengths, and is an expected experimental un- 
certainty. There are some similar features between our 



results and the experimental curve: there is a central re- 
gion sandwiched between two outer regions, like the low- 
field and the subcritical regions of the previous section. 
The size of the discontinuities grows towards lower fields 
at an apparently correct pace. The average period Ab 
for positive fields is longer than that for negative fields, 
and the ratio 1.30 : 1.17 seems to be of the right order. 
For fi — 60 and b ~ we predict a current of the order of 
0.6c$o[j,t/2(2irKR) 2 , where t is the thickness of the ring. 
Assuming t w 90nm and k ~ 0.5, this current amounts 
to 2 x 10 6 esu/sec. A current like this would induce a 
magnetic field of ~ 2G at the center of the ring; this is 
the correct order of magnitude. 

Figure [)] compares our results with a measurement in 
Ref. nil The micrograph at a corner shows an aluminum 
ring grown on a Hall cross which permits to measure tha 
average magnetic field produced by the supercurrent 
The dots are the experimental results. As can be seen 
from the micrograph, the width is close to 0.2i?, enabling 
a direct comparison with our results. The thickness is t rs 
0.12/im. The outer radius, as measured by optical means, 
is R = 1.1 ± 0.1/im. We picked R = 1.22/xm, in order 
to have the correct average periodicity. The temperature 
was 0.6±0.1K, which implies a coherence length between 
0.196/im and 0.217/xm and hence /i = 34 ± 3. The line 
in Fig. H was evaluated for /i = 30. The difference in 
height between the measured and the calculated results 
is not very significant, since several factors have been just 
estimated; however, the range of the central part is much 
wider than predicted. Again, the average period Ab for 
positive fields is longer than that for negative fields; this 
time the effect is more pronounced than in Ref. [l3] 

Since both experiments show a central part with two- 
fluxoid jumps which is considerably wider than in our 
results, we list some deviations from our idealizations 
that might give rise to this discrepancy. 

• We have neglected the induced field. The magnetic 
penetration depth A = nRyr^l 2 ~ 0.5 x 2300nm x 
60~ 1//2 Ri 150nm is only marginally larger than the 
thickness of the ring, so that our assumption is not 
clearly justified. The influence of the induced field 
on "superheating" for the case of thick disks has 
been considered in Ref. |7| and does not hint at a 
qualitatively different behavior. 

• An experimental ring never has perfect axial sym- 
metry. Thisjmperfection has been treated in previ- 
ous studiesr 8 !'! 39 ! and has a profound influence near 
the onset of superconductivity, but is not expected 
to be important deep inside the superconducting 
region. Large imperfections have .been found to in- 
hibit the size of the transitions.L3t!3 In the present 
problem, variation of the width with the height 
might also be important. 

• Decays most probably occur before the ideal 
metastability limit is reached, mainly due to elec- 
tromagnetic noise. When taking measurements, 
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noise can be filtered out by a lock-in amplifier, but ticle has benefited from comments raised by O. Megged, 
in order to avoid premature decays only shielding S. Pedersen and J. Rubinstein, 
could help. 
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